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Abstract. The young star velocity field is analysed by means of a galactic model which takes into account 
solar motion, differential galactic rotation and spiral arm kinematics. We use two samples of Hipparcos data, 
one containing O- and B-type stars and another one composed of Cepheid variable stars. The robustness of our 
method is tested through careful kinematic simulations. Our results show a galactic rotation curve with a classical 



value of A Oort constant for the O and B star sample (A 



13.7-13.8 km s" 1 kpc _1 ) and a higher value for 



Cepheids (A Cep = 14.9-16.9 km s _1 kpc _1 , depending on the cosmic distance scale chosen). The second-order 
term is found to be small, compatible with a zero value. The study of the residuals shows the need for a if -term 
up to a heliocentric distance of 4 kpc, obtaining a value K = —(1-3) km s _1 kpc~\ The results obtained for the 
spiral structure from O and B stars and Cepheids show good agreement. The Sun is located relatively near the 
minimum of the spiral perturbation potential (ipQ = 284-20°) and very near the corotation circle. The angular 
rotation velocity of the spiral pattern was found to be f2 p w 30 km s _1 kpc _1 . 

Key words. Galaxy: kinematics and dynamics - Galaxy: solar neighbourhood - Galaxy: structure - Stars: early- 
type - Stars: kinematics - Stars: variables: Cepheids 



• i— i 1. Introduction 

X 

$H 1 Young stars have been traditionally used as probes of the 
5^ i galactic structure in the solar neighbourhood. Their lu- 
minosity makes them visible at large distances from the 
Sun, and their age is not great compared to the dynamical 
evolution timescales of our galaxy. 

These stars show kinematic characteristics that can- 
not be explained by solar motion and differential galactic 
rotation alone. Small perturbations in the galactic grav- 
itational potential induce the formation of density waves 
that can explain part of the special kinematic features of 
young stars in the solar neighbourhood and also the ex- 
istence of spiral arms in our galaxy. The first complete 
mathematical formulation of this theory was done by Lin 
and his associates (Lin & Shu 1964 ; Lin ct al. 1969 ). 

Lin's theory has several free parameters that can be de- 
rived from observations. Two of the main parameters are 
the number of spiral arms (m) and their pitch angle (i). 
Although the original theory proposed a 2-armed spiral 



structure with i = —6°, as early as the mid- 70s Georgelin 
& Georgelin (1976) found 4 spiral arms with i — —12° 



from a study of the spatial distribution of HII regions. 
This controversy is still not resolved: in a review Vallee 
(1995) concluded that the most suitable value is m = 4, 
whereas Drimmel (2000) found that emission profiles of 
the galactic plane in the K band -which traces stellar 
emission- are consistent with a 2-armed pattern, whereas 
the 240 /im emission from dust is compatible with a 4- 



armed structure. In a recent paper, Lepine et al. (2001) 



described the spiral structure of our Galaxy in terms of a 
superposition of 2- and 4-armed wave harmonics, studying 
the kinematics of a sample of Cepheids stars and the l-v 
diagrams of HII regions. 

The angular rotation velocity of the spiral pattern £! p 
is another parameter of the galactic spiral structure. It 
determines the rotation velocity of the spiral structure as 
a rigid body. The classical value proposed by Lin et al. 



( 1969[ ) is fi p w 13.5 km s 1 kpc 1 . The angular rota- 
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tion velocity in the solar neighbourhood due to differen- 
tial galactic rotat ion is VLq s» 26 km s _1 kpc -1 (Kerr & 
Lynden-Bcll 1986 ). Thus, the value of fl p implies that the 
so-called corotation circle (the galactocentric radius where 
f2 p = f2) is in the outer region of our galaxy (G7 cor ~ 15-20 
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kpc, depending on the galactic rotation curve assumed). 
Nevertheless, several authors found higher values of O p , 
about 17-29 km s _1 kpc" 1 (Maroc hnik e t al. |l972j ; Creze 



& Mennessier [l973|; Byl fe Ovenden |1978 
Amaral & Lepinc 1997: Mishurov et al 
& Zenina 1999; Lepine et al 



Avedisova 



198E 



1997; Mishurov 



2001). These values place 



the Sun near the corotation circle, in a region where the 
difference between the galactic rotation velocity and the 
rotation of the spiral arms is small. This fact has very 
important consequences for the star formation rate in the 
solar neighbourhood, since the compression of the inter- 
stellar medium due to shock fronts induced by density 
wave s could be chiefly responsible for this process (Roberts 
1970|) . 



Other parameters of Lin's theory are the amplitudes 
of induced perturbation in the velocity (in the antigalac- 
tocentric and the galactic rotation directions) of the stars 
and gas, and the phase of the spiral structure at the Sun's 
position. The interarm distance and the phase of the spiral 
structure can be determined from optical and radio obser- 
vations (Burton |1971|; Bok & Bok |1974|; Schmidt-Kaler 



1975; Elmcgrccn 




The interarm distance gives us a 
relation between the number of arms and the pitch angle. 

In this paper we obtain the galactic kinematic param- 
eters from two samples of Hipparcos stars described in 
Sect. |^: one that contains O- and B-type stars, and an- 
other one composed of Cepheid variable stars. In Sect. || 
we propose a model of our galaxy, a generalization of that 
previously used by Comeron & Torra (1991). The authors 
only applied their model to radial velocities. The accu- 
rate astrometry of the Hipparcos satellite offers a good 
opportunity to also use the proper motion data. The res- 
olution of the condition equations, based on a weighted 
least squares fit, is explained in Sect. ||. An extensive set 
of simulations is performed in Sect. [s| in order to assess the 
capabilities of the method, that is, to analyze the influence 
of the observational errors and biases in the kinematic pa- 
rameters. We finally present our results and discussion in 
Sect. |6] and fj], respectively. 



2. The working samples 

2.1. Sample of and B stars 

Our initial sample contains 6922 Hipparcos O- and B-type 
stars. The astrometric data for these stars come from the 



Hipparcos Catalogue (ESA 1997), radial velocities from 
the compilation of Grenier (1997) and Stromgren pho- 
tometry from Hauck & Mermilliod's (1998) catalogue. 
The procedure followed to elaborate the sample consid- 



ered here is fully described in Torra et al. (2000; hereafter 
referred to as Paper I) along with a study of the possible 
biases in the trigonometric distances and the availabil- 
ity of radial velocities. Our sample contains 3915 stars 
with known distance and proper motions (2272 stars with 
known radial velocity). In Paper I we characterized the 
structure and kinematics of the Gould Belt system using 
this sample of O and B stars, establishing its boundary 
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Fig. 1. Star distribution in the X-Y galactic plane for the 
sample of O and B stars with 0.6 < R < 2 kpc. 

to a distance of about 0.6 kpc from the Sun. Taking into 
account the kinematic peculiarities of the Gould Belt, in 
this paper we do not consider those stars with R < 0.6 
kpc. In Fig. [l] we show the position of those stars with 
0.6 < R < 2 kpc projected on the galactic plane (X pos- 
itive towards the galactic center and Y positive towards 
the galactic rotation direction), which is our working sam- 
ple (448 stars; 307 of them with distance, radial velocity 
and proper motions and 141 with only distance and proper 
motions) . 



2.2. Sample of Cepheid stars 

The initial sample contains all the Hipparcos classi- 
cal Cepheids. Astrometric data were taken from the 



Hipparcos Catalogue (ESA 1997 ) , whe reas radial veloc- 
ities come from Pont et al. (|l994j , |l997| ) . 

Individual distances were computed following two 
period-luminosity (PL) relations. In both, periods come 
from the Hipparcos Catalogue (ESA 1997) and individual 
reddenings from Fernie et al.'s ( 1995[ ) compilation (contin- 
uous updating) . A classification between fundamental and 
overtone Cepheids from light curves and Fourier analysis 
was adopted (Beaulieu 1999] ) , using only the former in the 
kinematic analysis. 



The first PL relation (Luri 2000) adopts a slope from 
EROS (Beaulieu 1999) and corresponds to the short cos- 
mic distance scale: 



M 



Short 



= -1.08-2.721ogP 



(1) 



On the o ther h and, the second PL relation (Feast & 
Catchpole 1997 ) corresponds to the large distance scale: 



M Largo = _ L41 _ 3 4 g log p 



(2) 
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Fig. 2. Star distribution in the X-Y galactic plane for the 
sample of Cepheid stars with 0.6 < R < 4 kpc. Dis tance s 
computed from a short cosmic distance scale (Luri 200fj|) . 



The sample contains 186 stars with known distance and 
proper motions (165 stars with radial velocity). Their dis- 
tribution in the galactic plane (0.6 < R < 4 kpc) is shown 
in Fig. H (164 stars; 145 of them with distance, radial ve- 
locity and proper motions and 19 with only distance and 
proper motions). 

3. A galactic kinematic model for the solar 
neighbourhood 

3.1. Systematic velocity field in the galactic model 

We propose an extension of the bidimensional model pre- 
viously applied for radial velocities by Comeron & Torra 
( |l99l| ) . In our model we considered the systematic velocity 
components: 



v r = v Tl + v l2 + v l3 
V\ = v h + v h + v h 



(3) 



where subindex 1 refers to the solar motion contribution, 
subindex 2 to differential galactic rotation and subindex 3 
to spiral arm kinematics, respectively (see Appendix |Aj) . 
Solar motion is expressed through the three components 
of the Sun's velocity in galactic coordinates ({/©, Vq,Wq). 
Galactic rotation curve was developed up to second-order 
approximation, a r and b r being the first- and second-order 
terms, respectively: 



de 

dm 



e(tu Q ) 



2A 



dzu 2 



(4) 



where we show the relationship between a r and the A Oort 
constant. B Oort constant can be derived from: 



B = A 



<d{w 



Wq 



(5) 



We considered as free parameters the galactocentric dis- 
tance of the Sun (ujq) and the circular velocity at the 
Sun's position (0(n7©)). Finally, spiral arm kinematics 
was modelled within the framework of Lin's theory. We 
considered the number of the spiral arms (to) and their 
pitch angle (i) as free parameters, and we derived the 
perturbation velocity amplitudes in the antigalactoccntric 
and tangential directions (lib and 6b, respectively; they 
were considered as constant magnitudes, assuming that 
their variation with the galactocentric distance is smooth) , 
the phase of the spiral structure at the Sun's position (^0) 
and the parameter /©, which takes into account the dif- 
ference in the velocity dispersion between the solar-type 
stars and the considered stars. This parameter is defined 
by the relationship between the spiral perturbation veloc- 
ity amplitudes for the sample stars (lib, ©b) and for the 
Sun (n bQ , 6 b0 ): 



n b0 = \ " 2 + Xstais n h = / n b 



©b = 



1 - V 2 + Xq 

1 ^ "F ^stare 
1 — V 2 + Xr A 



-e h ee / Q e b 



(6) 



where v is the dimensionless rotation frequency of the 
spiral structure and x is the stability Toomre's number 



(Toomre 1969), which depends on the velocity dispersion 
of the considered stars (see details in Appendix |A]) . The 
inclusion of / Q is new with regard to the model proposed 
by Comeron & Torra ( |1991[ ). 
Eqs. (pi) can be expressed as: 



v T (R,l,b) = ^a,/j(i?,?,6) 
i=i 
10 

vi(R,l,b) = 5>,/](i?,Z,6) 



i=i 
10 



v h (R,l,b) = ^a J -/jW,&) 



(7) 



where the constants a,j contain combinations of the kine- 
matic parameters we wish to determine (Uq 7 V©, Wq, a T , 
b T , tpQ, lib, @b and /©) and fj(R,l,b) are functions of 
the heliocentric distance and the galactic longitude and 
latitude (see Eqs. (|A.21|) and (|A.22|)). 



3.2. Free parameters of our galactic model 

A 2-armed Galaxy was the first proposed view for our stel- 
lar system, mainly derived from HI and HII observations, 
but also from the spatial distribution of supergiant stars 
and other bright objects. These classical studies show the 
existence of at least two arms inside the solar circle (the 
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Sagittarius-Carina or —I arm and the Norma-Scutum or 
—II arm), one local arm (Orion-Cygnus or arm) and one 
external arm (Perseus or +1 arm). The Orion-Cygnus arm 
seems to be a local spur (Bok 1958). Lin et al. ( 1969] ) pro- 
posed a galactic system with 2 main spiral arms, where 
the Norma-Scutum and the Perseus arms are two seg- 
ments of the same arm. By taking into account the in- 
terarm distance between the Sagittarius-Carina and the 
Perseus arms, these authors deduced a pitch angle of —6°. 



But, as early as the mid-70s, Georgelin & Georgelin (1976) 
proposed a 4-armed galactic system with a pitch angle of 
— 12° from a study o f the spatial distribution of HII re- 
gions. However, Bash ( 1981 ) examined this 4-armed model 
and found that a 2-arm pattern predicts HII regions in the 
same direction and with the sam e radi al velocities as those 
used by Georgelin & Georgelin ( 1976 ) , provided that dis- 



persion velocities were considered. 

In some recent papers several authors have also called 
this classical view in question. Vallee ( 1995 ) reviewed the 
subject of the determination of the pitch angle and the 
number of spiral arms and concluded that the Galaxy has 
a pitch angle of i = — 12 ± 1° and that, taking into account 
the observed interarm distance, it would be a system of 4 
spiral arms. This is also the opinion expressed by Amaral 
& Lepine (1997), who, fitting the galactic rotation curve 
to a mass model of the Galaxy, found an autoconsistent 
solution with a system of 2 + 4 spiral arms (2 arms for 
2.8 < xn < 13 kpc and 4 arms for 6 < w < 11 kpc, with 
the Sun placed at w Q — 7.9 kpc) and a pitch angle of 
i = -14°. Englmaier & Gerhard ( |1999[ ) used the COBE 
NIR luminosity distribution and connected it with the 
kinematic observations of HI and molecular gas in l-v di- 
agrams. They found a 4-armed spiral pattern between the 
corota tion of the galactic bar and the solar circle. Drimmcl 
( 200C ) found that the galactic plane emission in the K 
band is consistent with a 2-armed structure, whereas the 
240 //m emission from dust is compati ble w ith a 4-armed 
pattern. In a recent work, Lepine et al. ( 2001 ) analyzed the 
kinematics of a sample of Cepheid stars and found the best 
fit for a model with a superposition of 2+ 4 spir al arms. 
Contrary to the model by Amaral & Lepine ( 1997 ), Lepine 
et al. allowed the phase of both spiral patterns to be inde- 
pendent, deriving pitch angles of approximately —6° and 
— 12° for m — 2 and m — 4, respectively. They argued that 
this spiral pattern is in good agreement with the l-v dia- 
grams obtained from observational HII data, though they 
admit that pure 2-armed model produces similar results. 
In the visible spiral structure of the Galaxy derived by 
Lepine et al. (see their Fig. 3) the Orion-Cygnus or local 
arm is seen as a major structure with a small pitch angle. 
In contrast, Olano (2001) proposed that the local arm is 
an elongated structure of only 4 kpc of length and a pitch 
angle of about —40° (in better agreement with the obser- 
vational determinations of the inclination of the local arm 
found in the literature) formed from a supercloud about 
100 Myr ago, when it entered into a major spiral arm. 

In the case of the galactocentric distance of the Sun 
and the circular velocity at the Sun's position there are 



also some inconsistencies among the different values found 
in the literature. In 1986, the IAU adopted the values 
vjq — 8.5 kpc and 9(-n7 Q ) = 220 km s _1 (Kerr & 
Lyndcn-Bcll 1986). Recently, several authors have found 
values of nearly 7.5 kpc for vjq (Racine & Harris 1989 
Maciel 1993). A complete review was done by Reid (1993) 



who concluded that the most suitable value seems to be 
vjq = 8.0 ± 0.5 kpc. In kinematic studies there are seri- 
ous discrepancies between different authors. Metzger et al. 
( |1998D found zu Q = 7.7±0 3 kpc and O(n7 ) = 237±12 km 
s _i , whereas Feast et al. ( 1998| ) found tuq = 8.5 ±0.3 kpc 
(both f rom r adial velocities of Cepheid stars). Glushkova 
et al. ( 1998 ) found Wq = 7.3 ± 0.3 kpc from a com- 
bined sample including open c luster s, red supergiants and 
Cepheids. Oiling & Merrificld ( 1998 ) calculated mass mod- 
els for the Galaxy and concluded that a consistent picture 
only emerges when considering tuq = 7.1 ± 0.4 kpc and 
0(n7 Q ) = 184 ± 8 km s _1 . This value for the galactocen- 
tric distance of the Sun is in very good agreement with 
the only direct distance determination (tuq = 7.2 ± 0.7 
kpc), which was made employing proper motions of H2O 
masers (Reid |l993[ ). 

In the view of all that, we decided to derive the kine- 
matic parameters of our model from different combina- 
tions of the free parameters involved. On the one hand, 
concerning spiral structure, two models of the Galaxy were 
considered: a first model with m = 2 and i — —6°, and a 
second one with m — 4 and i — —14°. Both models are 
consistent with an interarm distance of about 2.5-3 kpc, 
depending on the adopted value for the distance from the 
Sun to the galactic center. On the other hand, concerning 
the galactocentric distance of the Sun and the circular ve- 
locity at the Sun's position, two different cases were also 
taken into account: a first one with zuq = 8.5 kpc and 
Q(wq) = 220 km s , and a second one with zuq = 7.1 
kpc and O^©) = 184 km s _1 . In both cases, the angular 
rotation velocity at the Sun's position is Q@ — 25.9 km 
s _1 kpc -1 . 



4. Resolution of the condition equations 

We determined the kinematic parameters of the galactic 
model via least squares fit from the equations: 

10 



V* = 5> i £(JM,&) 
2=1 
10 

v\ = k R /ii cos b = ^ a,jfj(R, I, b) 

i=i 
10 

v h = kRpi h = 22ajfj(R,l,b) 



(8) 



where v r is the radial velocity of the star in km s" 1 , k 
= 4.741 km yr (s pc R is the heliocentric distance 

of the star in pc, /ii and /ib are the proper motion in 
galactic longitude and latitude of the star in " yr^ 1 and 
b the galactic latitude. To derive the parameters a,j, an 
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iterative scheme extensively explained in Fernandez (1998) 
was followed. 

The weight system was chosen as: 

1 



Pk 



k,ohs ^fc,cos 



(9) 



where <7 bs are the individual observational errors in each 
velocity component of the star, calculated by taking into 
account the correlations between the different variables 
provided by the Hipparcos Catalogue, and cr cos is the pro- 
jection of the cosmic velocity dispersion ellipsoid in the 
direction of the velocity component considered (see Paper 
I)- 

To check the quality of the least squares fits we con- 
sidered the x 2 statistics for N — M degrees of freeedom, 
defined as: 



A' 



fc=i 



[Vk - y(x k ;ai, ■ ■■,aio) 



(10) 



where x k are the independent data (sky coordinates and 
distances), y k the dependent data (radial and tangential 
velocity components), N the number of equations and M 
numbe r of p arameters to be fitted. According to Press 
et al. ( 1992| ), if the uncertainties (cosmic dispersion and 
observational errors) are well-estimated, the value of x 2 
for a moderately good fit would be \ 2 ~ N ~ with an 
uncertainty of yj2 (N - M). 

As we did in Paper I, to eliminate the possible 
outliers present in the sample due to both the exis- 



tence of high residual velocity stars (Royer 1997) or 
stars with unknown large observational errors, we re- 
jected those equations with a residual larger than 3 times 
the root mean square residual of the fit (computed as 
\/[yk — y(xk', ai, o,iq)] 2 /N) and recomputed a new set 
of parameters. 

5. Test of robustness 

The spiral arm potential is expected to contribute between 
5-10% to the whole galactic gravitational field in the so- 
lar neighbourhood. This small contribution, and the large 
observational errors and constraints present on the spatial 
and kinematic parameters of distant stars, make it very 
difficult to quantify the kinematic perturbations induced 
by the spiral potential. Then, the results found in the liter- 
ature have been characterised by significant uncertainties 
and discrepancies. A good example of this is the contradic- 
tory results obtained for the phase of the spiral structure 
at the Sun's position (see Sect. 7.2): between arms or near 
an arm? Interesting questions to answer after the release 
of the Hipparcos data are: 

— How does the unprecedented astrometric precision pro- 
vided by Hipparcos help to diminish such uncertain- 
ties? Are they small enough to measure the spiral arm 
kinematic parameters to any useful degree of accuracy? 

— Can we quantify the biases induced by our observa- 
tional errors and constraints? 



Table 1. Number of stars with distance and proper mo- 
tions (in brackets those stars with also radial velocity) 
in several distance intervals. Distances for Cepheids com- 
puted from Luri's (200C) PL relation. 



sample of O and B stars 



0.1 < R < 2 kpc 0.6 < R < 2 kpc 



3418 (1903) 



448 (307) 



sample of Cepheid stars 



0.1 < R < 4 kpc 0.6 < R < 2 kpc 0.6 < R < 4 kpc 



119 (111) 



103 (95) 



164 (145) 



— How can the correlations among variables present in 
the condition equations affect the kinematic parame- 
ters? 

With regard to our stellar data, as seen in Sect. 0, both 
O and B star and Cepheid samples suffer from different 
observational limitations. Although the O and B star sam- 
ple is large in number, it is limited in distance to no more 
than 1.5-2 kpc from the Sun (see Fig. [l] and Table |l|). In 
contrast, the Cepheid sample reaches distances of up to 
about 4 kpc, but the number of stars with reliable data 
still remains very small (Fig. | and Table Q). 

Numerical simulations were performed to answer the 
above questions, that is, to evaluate and quantify all the 
uncertainties and biases involved in our resolution pro- 
cess. In Appendix [b| we present the detailed procedure 
followed to build simulated samples as similar as possible 
to the real data for both O and B stars and Cepheids. 
This Appendix also contains an exhaustive analysis of the 
full set of cases which were simulated. Next we present the 
main conclusions arising from this work, which substan- 
tially contribute to the analysis of the real data: 

— As a general trend, owing to the correlations between 
variables, the biases and uncertainties in some kine- 
matic parameters substantially vary when changing 
the real values of i/jq and £! p . 

— The first- and second-order terms of the galactic rota- 
tion curve are accurately derived, with an uncertainty 
of less than 1.3 km s _1 kpc~ x (km s _1 kpc -2 ) in all 
cases. Whereas a negative bias for both parameters is 
detected for O and B stars (up to —0.7 km s _1 kpc -1 
for a T and —1 km s _1 kpc -2 for b T ), it is negligible for 
Cepheids. 

— For the spiral structure parameters, we obtain better 
results from O and B stars than from Cepheids. This 
is due to the larger number of O and B stars. For both 
samples we found a clear dependence of some param- 
eters with the assumed values for ipQ and f2 p . We can 
expect a bias in ipo up to ±20°, and uncertainties of 
10-20° for O and B stars and 30-60° for Cepheids. For 
Hb and 0b we found biases up to ±2 km s _1 and uncer- 
tainties of 1-2 km s~ x . We did not find an important 
bias for fi p neither for O and B stars nor Cepheids, 
though its standard deviation is large (4-10 km s _1 
kpc -1 for O and B stars and 10-20 km s _1 kpc -1 for 
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Cepheids). These ranges of values have allowed us to 
check how compatible are the results independently 
obtained for O and B stars and Cepheids. 
— If we choose an incorrect set of free parameters (m, 
i, tuq, Q(vuq)), results change, but not to a large ex- 
tent. This change is similar to the dispersion obtained 
when solving the condition equations from the differ- 
ent simulated samples (see Appendix ^ for details). 
Therefore, it will be difficult to decide between 2- and 
4-armed models. 

From these results, in Appendix [b| we conclude that, 
with the present available observational data, we are able 
to determine the kinematic parameters of the galactic 
model proposed in this paper, though we are not able to 
decide between a 2- or 4-armed Galaxy. 

6. Results 

In this section we present the results obtained from our 
working samples. To avoid those stars belonging to the 
Gould Belt, which can produce important deviations in 
our results (see Paper I), we always considered stars fur- 
ther than 0.6 kpc from the Sun. Following the results ob- 
tained in Appendix |b[ the working distance intervals are 
0.6 < R < 2 kpc for O and B stars and 0.6 < R < 4 kpc 
for Cepheids. 

A first set of results is presented in Tabl e c onsider- 
ing a classic al view of our galaxy (Lin et al. 1969t Kerr & 
Lynden-Bcll 1986 ): we supposed a differential galactic ro- 
tation (a r , 6 r ) and a system with 2 spiral arms and a pitch 
angle of —6°, the Sun placed at a galactocentric distance 
of 8.5 kpc a the circular velocity at the Sun's position of 
220 km s _1 . 

Although we do not only present solutions from radial 
velocity or proper motion data, we would like to remark 
that the parameters obtained when solving these cases 
(and comparing with the combined solutions presented in 
Table ||) are compatible between themselves within the 
errors bars. This is true for the sample of O and B stars 
and also in the case of the sample of Cepheids (for both 
short and large cosmic scales). Nevertheless, as we deter- 
mined in Paper I, for this kind of kinematic analysis it is 
advisable to solve a combined resolution, in order to mini- 
mize the influence of the correlations present between the 
different parameters. This is especially important in the 
present case, since some correlations reach a large value 
(up to 0.8, though in the majority of cases they do not 
exceed 0.3). In Appendix [b| we check that these correla- 
tions do not impede our obtaining reliable results from 
our stellar samples. 

If we compare the results obtained from O and B stars 
and Cepheids, we observe several discrepancies: differences 
up to 3.5 km s _1 kpc -1 in the A Oort constant (derived 
from a r ) and up to 125° for the phase of the spiral struc- 
ture at the Sun's position. We must keep in mind that, 
whereas this parameter (like f2 p and w COT ) is independent 
of the sample used, lib, ©b and Jq depend on the cosmic 



Table 2. Resolution for the samples of O and B stars 
and Cepheids (SCS: short cosmic scale; LCS: large cos- 
mic scale). Units: Uq, Vq, Wq, li^, Ob and a in km 
s -1 ; a r in km s _1 kpc -1 ; b r in km s _1 kpc -2 ; ip Q in 
degrees. x 2 /N is the value of % 2 divided by the num- 
ber of equations minus the degrees of freedom. A cos- 
mic dispersion of {ay, &w) = (8,8,5) km s _1 and 
(crj/, ay , aw) = (13,13,6) km s _1 was used for O and 
B stars and Cepheids, respectively. 





O and B stars 


Cepheid stars 




0.6 < R < 2 kpc 


0.6 < R < 4 kpc 






SCS distances 


LCS distances 


u e 


8.8 ± 0.7 


6.5 ± 1.2 


8.3 ± 1.2 


Vq 


12.4 ± 1.0 


10.4 ± 1.9 


9.3 ± 2.1 


Wq 


8.4 ± 0.5 


5.7 ± 0.7 


7.0 ± 0.8 


GSr 


-1.3 ± 1.0 


-8.1 ± 1.2 


-4.5 ± 1.1 


6r 


-0.8 ± 1.5 


1.8 ± 0.8 


-0.1 ± 0.8 


IpQ 


45. ± 52. 


282. ± 20. 


306. ± 24. 


n b 


1.6 ± 1.1 


0.1 ± 1.6 


0.7 ± 1.4 


e„ 


2.1 ± 1.5 


4.9 ± 1.7 


4.8 ± 1.7 


/© 


0.42 ± 0.10 


0.96 ± 0.01 


0.96 ± 0.01 




12.11 


11.57 


11.77 




2.07 


0.84 


0.84 



dispersion of the sample (thus, they do not have the same 
value for O and B stars and Cepheids). 

As in Paper I, we found values of x 2 /N ~ 2 for O and 
B stars. We think that the difference from the expected 
value (x 2 /N « 1) could be due to an underestimation 
of the errors in the photometric distances and/or radial 
velocities for far stars. In the case of Cepheids we found 
values which agree with the expected x 2 /N ~ 1. 

When studying the residuals of the equations we re- 
alized that radial velocity equations have a non-null av- 
erage residual of about —(3-4) km s" 1 for both O and 
B stars and Cepheids. In the case of O and B stars (see 
Fig. ||) we can observe the peculiar motion of some OB 
associations. The kinematics of those associations located 
near the Sagittarius arm (e.g. Sgr OBI, Ser OBI and Set 
OB2, placed at 0° < I < 30°) were studied by Mel'nik 
et al. ( 1998 ), who found that their motions are in gen- 
eral agreement with what would be expected according to 
Lin's theory if these stars are located within the corota- 
tion radius. But we found these stars have a large residual, 
even taking into account in our model the galactic spiral 
structure, possibly due to the fact that they are still re- 
flecting a motion peculiar to their birthplaces owing to 
their youth. Another region with a high residual is that 
located in the direction of the open clusters h and x P er 
(NGC 869 and 884), at I » 135°, where we observe a 
group of stars (distributed in an area of 10° in the sky) 
with a high residual directed towards the Sun's position. 
Thus, our results indicate that these stars do not fit the 
systematic velocity field defined by the whole sample. In 
this region of the sky the associations Per OBI and Cas 
OB6 are to be found. The former is generally considered 
to include h and x P er an( i surrounding O and B stars. 
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However, the stars in Fig. ^ have distances 1.1 < R < 1.4 
kpc, whereas the open clusters are located further on, at 
R pa 2-2.5 kpc (Schild 1967). In our initial sample (see 
Sect. ^|), we have in this region a group of stars with dis- 
tances 2.0 < R < 2.6 kpc, most of them belonging to Per 
OBI according to the membership list done by Garmany & 



Stencel (1992). But we did not use these stars in our kine- 
matic analysis, due to their large distance errors. Garmany 
& Stencel (1992) questioned if h and x P er really belong 
to Per OBI due to the fact that B main sequence stars 
do not define a ZAMS at the distance of the open clus- 
ters. These authors derived a distance modulus of 11.8 
(R = 2.3 kpc) for Per OBI and 11.9 (R = 2.4 kp c) for 
Cas OB6. More recently, Mel'nik & Efremov fll995[ ) stud- 
ied the spatial distribution of O and B stars within 3 kpc 
from the Sun and derive d a ne w partition into OB associ- 
ations using Battinelli's ( 1991 ) modification of the cluster 
analysis method. They found that Per OBI split into four 
groups, two at R « 1.6 kpc and another two at R rj 1.9 
kpc. Two of them are reliable at a 90% confidence level 
(one at 1.6 kpc and another at 1.9 kpc), with a mean he- 
liocentric radial velocity of —21.6 km s _1 and —41.4 km 
s^ 1 , respectively. The first association found by Mel'nik & 
Efremov might correspond to the group of stars we have 
at 1.1 < R < 1.4 kpc. In view of all that, the physical 
explanation of the large negative radial residuals present 
in this region, even considering the contribution from the 
spiral arm kinematics, is still unresolved. 

If we exclude those stars belonging to both regions 
from our calculations, the average residual in radial ve- 
locity decreases to —(1-2) km s _1 . For Cepheids, as men- 
tioned above, the residuals have the same trend and reach 
up to —(3-4) km s -1 , though in this case we cannot clearly 
identify groups of stars with a common motion. Therefore, 
the mean negative residual motion observed does not seem 
to be explained as the peculiar motion of some groups 
of stars. Moreover, we found that, for both samples, this 
residual motion apparently does not depend on the galatic 
longitude of the stars. 

Taking into account that Cepheids, with a large mean 
distance, show larger residuals than O and B stars, we 
can suppose that this systematic trend could be explained 
through the addition of a K-tevm in Eqs. (^), as developed 
in Paper I. Then, the following terms can be added to Eqs. 



fix 
fix 
f h xx 



R cos 2 b 


— i?sinocosf> 



(11) 



which allow us to derive the parameter an = K. Table 
|3| shows the results obtained in this way, and considering 
different sets of free parameters (as proposed in Sect. |3.2| ). 
The correlation coefficient between K and the other pa- 
rameters is always low, being smaller than 0.4 for O and 
B stars and 0.2 for Cepheids. As stated in Appendix [bJ a 
change in the value of the free parameters (m, i, vuq and 
Q(v&q)) does not significantly alter the kinematic param- 



o - 



GR - 



- 7 v / t 



GC 



1 



\ 



= 25 km/s 



Y (kpc) 

Fig. 3. Residual heliocentric space velocity vectors pro- 
jected on the galactic plane for O and B stars. GR stands 
for galactic rotation direction and GC for galactic center 
direction. 

eters derived. In Table || we also show the results when 
rejecting those stars belonging to the problematic regions 
mentioned above in the case of O and B stars. We observe 
that there is not a great difference in the results between 
these cases, which are always compatible within the er- 
rors. So, we conclude that the observed peculiar motion 
of some OB associations present in our sample does not 
alter the main parameters of our model. In view of that, 
in the next section we discuss the results obtained when 
considering all the stars, without rejecting associations. 

7. Discussion 

7.1. Galactic rotation curve 

A difference in the A Oort constant of approximately 3 
km s _1 kpc -1 between solutions for O and B stars and 
Cepheids with short cosmic scale distances is obtained, 
whereas large cosmic scale provides an intermediate value: 



^OB 

A Ccp 
^Short 



13.7-13.8 km s _1 kpc -1 
16.6-16.9 km s _1 kpc -1 



(12) 



A Sr P g o ~ 14-9-15.1 km s" 1 kpc" 1 

As can be seen in Appendix |b], these differences cannot be 
explained by the observational constraints present in both 
samples. A similar discrepancy was found by Frink et al. 
( P-996D , who derived A = 14.0 ± 1.2 km s _1 kpc -1 from a 
sample of O and B stars, and A — 15.8±1.6 km s _1 kpc -1 
from a sample of Cepheids (in both cases the authors only 
considered those stars with a heliocentric distance of less 
than 1 kpc). Nonetheless, our values for Cepheids d o not 
reach the higher values obtained by Glushkova et al. ( |l998| ; 
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Table 3. Resolution by least squares fit considering differ- 
ent values of the imposed galactic parameters m, i, Wq and 
Q(wq). Case A: m = 2, i = — 6°, w & = 8.5 kpc, 6(^7©) = 
220 km s -1 ; Case D: m = 4,i = ~U°,w Q = 7.1 kpc, 
0(tz7 ) = 184 km s -1 . Units: U e , U , W e , II b , O b and 
a in km s ; a T , K, A, B and f2 p in km s -1 kpc -1 ; b T in 
km s -1 kpc -2 ; xp Q in degrees; Aw COI = w Q — w cot in kpc. 
X 2 /N is the value of \ 2 divided by the number of equa- 
tions minus the degrees of freedom. Gl and G2 stand for 
those groups of stars located in < I < 50°, 1 < R < 2 
kpc (29 stars) and 130 < Z < 140°, 1 < R < 2 kpc (14 
stars), respectively. 

O and B stars with 0.6 < R < 2 kpc 





Case A 


Case D 


Case D 


Case D 


excluding Gl excluding G2 


Uq 


9.2 ± 0.7 


10.0 ± 0.7 


9.2 ± 0.8 


9.6 ± 0.8 


V e 


12.7 ± 1.1 


12.7 ± 1.0 


13.2 ± 1.0 


12.3 ± 1.1 


Wq 


8.3 ± 0.5 


8.3 ± 0.5 


8.3 ± 0.5 


8.4 ± 0.5 


a r 


-1.7 ± 1.0 


-1.5 ± 1.0 


-2.3 ± 1.2 


-0.9 ± 1.0 


b r 


-0.4 ± 1.5 


0.5 ± 1.3 


2.3 ± 1.5 


0.0 ± 1.4 


K 


-3.2 ± 0.7 


-2.8 ± 0.7 


-1.4 ± 0.7 


-2.4 ± 0.7 


V'Q 


20. ± 53. 


329. ± 47. 


315. ± 32. 


8. ± 52. 


n b 


3.1 ± 1.1 


2.6 ± 1.0 


3.2 ± 1.1 


2.8 ± 1.1 


e b 


2.0 ± 1.3 


2.1 ± 1.3 


3.2 ± 1.5 


1.8 ± 1.3 


/© 


0.23 ± 0.57 


0.35 ± 0.25 


0.40 ± 0.10 


0.19 ± 0.77 


A 


13.8 ± 0.5 


13.7 ± 0.5 


14.1 ± 0.6 


13.4 ± 0.5 


B 


-12.7 ± 0.8 


-13.2 ± 0.7 


-14.1 ± 0.8 


-13.0 ± 0.7 


ftp 


45. ± 14. 


33. ± 6. 


32. ± 3. 


36. ± 8. 


Azu COT 


3.6 ± 1.6 


1.5 ± 0.9 


1.1 ± 0.6 


2.0 ± 1.2 




11.85 


11.88 


11.74 


11.85 


xVn 


1.89 


1.90 


1.90 


1.89 


Cepheid stars with 0.6 < R < 4 kpc 




SCS distances 


LCS distances 




Case A 


Case D 


Case A 


Case D 


Uq 


6.5 ± 1.2 


6.5 ± 1.2 


8.4 ± 1.2 


8.2 ± 1.2 


Vq 


10.5 ± 1.9 


11.0 ± 2.1 


9.1 ± 2.1 


10.1 ± 2.3 


Wq 


5.7 ± 0.7 


5.7 ± 0.8 


7.0 ± 0.8 


7.0 ± 0.8 


a r 


-7.9 ± 1.2 


-7.4 ± 1.2 


-4.4 ± 1.1 


-3.9 ± 1.1 


b t 


1.6 ± 0.8 


0.5 ± 0.8 


0.1 ± 0.8 


-0.9 ± 0.8 


K 


-0.8 ± 0.5 


-1.0 ± 0.5 


-1.2 ± 0.5 


-1.2 ± 0.5 




286. ± 21. 


284. ± 39. 


310. ± 21. 


321. ± 43. 


n b 


0.0 ± 1.6 


-1.4 ± 1.6 


0.9 ± 1.3 


0.4 ± 1.2 


e b 


4.8 ± 1.7 


2.4 ± 1.7 


5.6 ± 1.7 


2.7 ± 1.7 


fe 


0.96 ± 0.01 


0.95 ± 0.01 


0.96 ± 0.01 


0.96 ± 0.01 


A 


16.9 ± 0.6 


16.6 ± 0.6 


15.1 ± 0.6 


14.9 ± 0.6 


B 


-13.7 ± 0.4 


-13.2 ± 0.4 


-13.0 ± 0.4 


-12.5 ± 0.4 


ftp 


26. ± 3. 


23. ± 4. 


28. ± 3. 


27. ± 2. 


AzUcor 


0.0 ± 1.0 


-0.7 ± 1.2 


0.5 ± 0.8 


0.2 ± 0.6 




11.56 


11.66 


11.70 


11.90 




0.84 


0.85 


0.83 


0.86 



kpc 1 ), Mishurov et al. ( 
kpc -1 ), Mishurov & Zcnina 



1997 



1999 



A = 19.5 ± 0.5 km 
A = 20.9 ± 1.2 km s" 

A = 18.8±1.3kms -1 kpc -1 ) and Lcpincct al. (|200l[ ; A = 
17.5 ± 0.8 k m s -1 kpc -1 ). Pont et al. fll994j ) and Metzger 
et al. (199S), from radial velocities of Cepheid stars, found 
values of A = 15.9 ± 0.3 km s -1 kpc -2 and A = 15.5 ± 
0.4 km s -1 kpc -2 , respectively. More recently, and using 



proper motions and distance calibration from Hipparcos 
data on Cepheid stars, Feast & Whitelock (1997) found 
a value of A = 14.8 ± 0.8 km s -1 kpc -1 . Using a similar 
sample, Feast et al. ( |1998| ) found A = 15.1 ± 0.3 km s -1 
kpc -1 from radial velocities. As we have mentioned in 
Sect. we found good coherence for radial velocity, proper 
motion and combined solutions for both cosmic distance 
scales. 

An attempt to explain these discrepancies was made 
by Oiling & Merrifield ( 1998D , who studied the variation 
of the A and B Oort junctions and found that they sig- 
nificantly differ from the general ~ 0(w)/w dependence 
expected for a nearly flat rotation curve. Inside the so- 
lar circle, the value of A rises to 18 km s -1 kpc -1 for 
Aw = w — vjq ps —0.5 kpc, disminishes to 16 km s -1 
kpc -1 for Aw w —1.2 kpc, and rises continuously for 
Aw —1.5 kpc, to 19 km s -1 kpc -1 for Aw w —2 kpc 
(see Fig. 3 in Oiling & Merrifield |1998[ ). Contrary to that, 
beyond the solar circle A decreases to 10-12 km s -1 kpc -1 , 
maintaining this value in the interval Ss Aw <, 2.5 kpc. 
Although A is a local parameter, describing the local shape 
of the rotation curve, Oiling & Merrifield already pointed 
out that the discrepancies in the results published in the 
literature may be produced by their dependence on the 
galactocentric distance. Our O and B stars are distributed 
along all the galactic longitudes, whereas the Cepheids are 
predominantly concentrated inside the solar circle, with 
a peak in the spatial distribution for Aw k, —0.6 kpc 
corresponding to the Sagittarius-Carina arm (see Figs, [l] 
and |2|). For O and B stars we found a classical value of 
A w 14 km s -1 kpc -1 (a value between 10-12 and 16-18 
km s -1 kpc -1 ), whereas for Cepheids a value A k, 15-17 
km s -1 kpc -1 (depending on the PL relation considered) 
was derived, in agreement with Oiling & Merrifield's as- 
sumptions. 

From the results obtained in Appendix^, we would ex- 
pect uncertainties in the second-order term of the rotation 
curve of « 1.0 km s -1 kpc -2 for O and B stars and ~ 0.5 
km s -1 kpc -2 in the case of Cepheids. Taking this and the 
results in Table || into account, we can state that b r does 
not differ from a null value more than 2 km s -1 kpc -2 . 
Pont et al. Q1994D found a v alue b T = -1.7 ± 0.2 km s -1 
kpc -2 , whereas Feast et al. (1998) found b T = —1.6 ± 0.2 
km s -1 kpc -2 , both using radial velocities of Cepheid stars 
(the latter with a Hipparcos distance calibration) . A large 
positive value was found by Lepine et al. ( |200l| ), who de- 
rived b r = 5.0 ± 1.0 km s -1 kpc -2 from their sample of 
Cepheid stars. As we will see in Sect. [7.2| , these authors 
also found a large value for the amplitude of the velocity 
component in the galactic rotation direction due to the 
spiral potential (O b ). Without specific simulations, con- 
sidering both their observational data and resolution pro- 
cedure, it is difficult to guess how the correlations between 
both parameters can affect their determination. 
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Fig. 4. Star distribution in the X-Y galactic plane for O 
and B stars (filled circles) with 0.6 < R < 2 kpc and 
Cepheids (empty circles) with 0.6 < R < 4 kpc (short 
cosmic scale distances). Spiral arms were drawn consider- 
ing ipQ = 330°. Dotted lines show the center of the spiral 
arms (ip = 0°) and dashed lines draw their approximate 
edges (tp — ±90°). Looking at this figure we must have in 
mind the possible drift between the position of the optical 
tracers in the spiral arms (i.e., our young stars) and the 
minimum of the spiral potential (ip = 0°). 



7.2. Spiral structure 

Appendix [b| demonstrates that the available observational 
data allow the characterization of the galactic spiral struc- 
ture, though the biases and uncertainties on the parame- 
ters have to be taken into account in the interpretation of 
the results. Furthermore, we realized that it is very dim- 
cult to establish the correct number of spiral arms of the 
Galaxy. 

A first remarkable result is the fairly good coherence 
obtained for the phase of the spiral structure at the Sun's 
position ipQ when using different free parameters (cases 
A, D) or different samples compared to the great discrep- 
ancies found in the literature (see below). We take into 
consideration that lib, Ob and depend on the cosmic 
dispersion and so they do not have the same value for O 
and B stars and Cepheids. 

From Fig. and assuming that Cepheids fairly trace 
the center of the Sagittarius-Carina arm, we can infer a 
value of ipo w 250° (the center of the inner visible spiral 
arm -the Sagittarius-Carina arm- at about 1 kpc from 
the Sun), depending on the exact value of the interarm 
distance. Nevertheless, the density wave theory predicts 
that the center of the visible arm (traced by its young 
stars) does not coincide wit h the posi tion of the spiral 
potential minimum (Roberts ^969] |1970| ). We must bear in 



mind that the kinematically derived value for ipQ informs 
us about the position of the Sun with respect this potential 
minimum, not with respect the visible arm. In Table ^ we 
found values inside the interval: 



284-20 c 



(13) 



We take into consideration that in the minimum of the spi- 
ral potential (near the center of an arm) -0 = 0°, whereas 
in the inner edge of the arm tp w 90° and in the outer 
edge ijj w —90° = 270°. In the case of the phase of the 
spiral structure at the Sun's position we expect a bias 
of AV>© = V>o btainod - 0Q al ~ 20° for O and B stars 
and AipQ w 0° for Cepheids, with uncertainties of « 25° 
and w 75°, respectively (see Appendix ^). These high 
uncertainties can explain the range of values obtained. 
According to the value found for ipQ, the Sun is located 
between the center and the outer edge of an arm, nearer 
to the former (see Fig. |J) . 

Our result stands in apparent contradiction to the clas- 
sical mapping of the spiral structure tracers (Schmidt- 
Kaler 1975 Elmegreen [1985 ), which locates the Sun in a 
middle position between the inner arm (Sagittarius-Carina 
arm) and the outer one (Perseus arm), that is, -0 Q ss 180°. 
The local arm (Orion-Cygnus arm) is normally considered 
as a local spur rather than a real arm. In our model, as 
we supposed an interarm distance of approximately 3 kpc, 
the local arm is also considered a local spur (see Sect. 3.2). 
But the value we obtained for ipo differs significantly from 
180°. If we consider a difference of 30-100° between the op- 
tical tracers and the potential minimum of the spiral arm, 
we notice that our ipQ value is in good agreement with 
the picture of the spiral structure that emerges from the 
spatial distribution of stars in Fig. ^, with the inner spi- 
ral arm at approximately 1 kpc from t he Su n. A similar 
result was obtained by Mel'nik et al. ( 1998| ), who found 
that nearly 70% of the stars in the OB associations of the 
Sagittarius-Carina arm have a residual motion (after cor- 
recting their heliocentric velocities for solar motion and 
galactic rotation) in the direction opposite to the galac- 
tic rotation, as one would expect for those stars between 
the inner edge and the center of the arm. Then, they also 
found a shift between the optical position of the visible 
arm (traced by young stars) and the minimum of the spi- 
ral potential. Therefore, from our results we conclude that 
the Sun is placed relatively near the potential minimum 
of the Sagittarius-Carina arm, and the Perseus arm is lo- 
cated far away, at about 2.5 kpc from the Sun. 

Our range of values for ipQ includes that obtained 



by Creze & Mennessier (1973) from a sample of 0-B3 
stars, ipo — 352 ± 30°. Later, Mennessier & Creze (1975) 
found from a sample of 0-B3 stars a value of ipo « 90°, 
which locates the Sun near the inner edge of an arm. 
Other autho rs ob tained contradictory results. Gomez & 
Mennessier ( 1977 ) found the Sun's position near the edge 
of an arm from several samples of stars from FK4 and 
FK4 Supplement Catalo gues. Byl & Ovenden (|l978j) and 
Comeron & Torra ( 1991 ) obtained values of tp® = 165±1° 
and -00 = 135 ± 18° respectively, from samples of O and 
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B stars. Mishurov et al. ( |1997| ) found ip Q = 290 ± 16° 
from a sample of Cepheid stars with radial velocities. 
Mishurov & Zenina ( |1999[ ) found ip & = 322 ± 9° (sup- 
posing m — 2 and Wq = 7.5 kpc), from the same sample 
of Cepheid stars, but also including the Hipparcos proper 
motions. When these authors supposed m = 4 they found 
^ Q = 340 ± 9°. More recently, Rastorguev et al. fl2001| ) 
found ipQ — 274 ± 22° from a sample of 55 open clus- 
ters younger than 40 Myr and 67 Cepheids with periods 
smaller than 9 days, all of them at R < 4 kpc. We can see 
that there is poor agreement between the results published 
in the literature, though the last ones (the only ones with 
Hipparcos data) are all included in our range of possible 
values for ipQ. 

The velocity amplitudes due to the spiral perturbation 
obtained are less than 4 km s" 1 for O and B stars (lib ~ 3 
km s _1 , 0b ~ 2-3 km s" 1 ) and 6 km s" 1 for Cepheid 
stars (lib w -1-1 km s" 1 , 9 b « 2-6 km s" 1 ). The dif- 
ference between results for O and B stars and Cepheids 
is explained by Lin's theory as a result of their slight de- 
pendence on the galactocentric distance and the differ- 
ent cosmic dispersion for both samples. Mishurov et al. 
( |1997| ) found n b = 6.3 ± 2.4 km s" 1 and 6 b = 4.4 ± 2.4 
km s _1 from their radial velocity data of Cepheid stars 
within 4 kpc from the Sun, supposing a 2-armed spiral 
pattern. Mel'nik et al. fll999| ) found n b = 6.4 ± 1.2 km 
s" 1 and Ob = 2.4 ± 1.2 km s" 1 from Cepheids within 
3 kpc from the Sun. Mishurov & Zenina (1999) found 
n b = 3.3 ± 1.6 km s" 1 and 6 b = 7.9 ± 2.0 km s" 1 for 
m = 2 and n b = 3.5 ± 1.7 km s" 1 and 9 b = 7.5 ± 1 8 km 
s" 1 for m = 4. From the same sample, Lepine et al. ( |2001 ) 
found n^ 2 = 0.4 ± 3.0 km s" 1 , Q r £= 2 = 14.0 ± 3.0 km 



with the classical fi p « 13.5 km s 1 kpc 1 proposed by 



s" 1 and n™= 4 = 0.8 ± 3.3 km s 



l Q m=4 = 10 .9±2.9 km 



s . Their 2+4-armed model yields large values for 9b, 
implying a large value for the ratio between the spiral po- 
tential and the axisymmetric galactic field, much greater 
than that 5-10% normally accepted. According to the re- 
sults obtained in Appendix |b], the present observational 
uncertainties, biases and the correlations involved in the 
resolution procedure cannot completely explain the dis- 
crepancies among the values found in the literature. They 
can be attributed to our poor knowledge of the real wave 
harmonic structure of the galactic spiral pattern or to the 
approximation performed in the linear density-wave the- 
ory. 

The derived angular rotation velocity of the spiral pat- 
tern is: 



flp ~ 30 km s 1 kpc 1 



(14) 



though there is a great dispersion around this value in our 
results from different samples and cases (dispersion of 2-7 
km s" 1 kpc" 1 , but up to 15 km s" 1 kpc -1 in one extreme 
case). This is an expected dispersion, taking into account 
the results obtained in Appendix [b| (standard deviation of 
2-5 km s" 1 kpc -1 for O and B stars, but up to 15 km s _1 
kpc -1 for ipQ k 270°; for Cepheids dispersions are higher, 
of about 10-20 km s" 1 kpc -1 ). Such a value places the 
Sun very near the corotation circle and is not in agreement 



Lin et al. (1969). Nevertheless, other studies show results 
similar to ours. Avedisova (198E) obtained £l p = 26.8 ± 
2 km s" 1 kpc -1 from the spatial distribution of young 
objects of different ages in the Sagittarius-Carina arm. 
More recently, Amaral & Lepine (1997), working with a 
selection of young members of the open cluster catalogue 



by Mermilliod (1986), obtained a value of Q p w 20-22 km 



s 1 kpc 1 . Mishurov et al. (1997) and Mishurov & Zenina 
( 1999| ) used their sample of Cepheid stars and found fl p = 
28.1 ± 2.0 km s" 1 kpc -1 andf> p « 27.7 km s" 1 kpc" 1 , 
respectively. Lepine et al. (|200l|) found Q p n = 2 « ft™= 4 » 
26.5 km s _1 kpc" 1 from their sample of Cepheid stars. 



Rastorguev et al. ( 2001 ) thought that the evaluation of 



fl p from kinematical data alone cannot be resolved. But 
our simulations (see Appendix |b|) seem to indicate that, at 
least, the tendency to find high values of fi p is confirmed, 
though there is still an uncertainty of about 5-10 km s" 1 
kpc -1 in its value. These large values for Q p can explain in 
a natural way the presence of a gap in the galactic gaseous 
disk (see Kerr 1969 and Burton 1976 for observational 



evidences and Lepine et al. 2001 for simulated results), 
since they placed the Sun near the corotation circle, where 
the gas is pumped out under the influence of the spiral 
potential. 



7.3. The K-term 

In relation to the K-term, we found good agreement be- 
tween both O and B stars and Cepheids. A value of 
K w —(1-3) km s _1 kpc -1 is found in all cases, confirm- 
ing an apparent compression of the solar neighbourhood 
of up to 3-4 kpc. In our opinion, what most clearly proves 
the existence of a non-null value of K is that it is inde- 
pendently found from the samples of O and B stars and 
Cepheids, which have a very different spatial distribution, 
an independent distance estimation and a different way of 
deriving their radial velocities. As we can see comparing 
Tables || and || the inclusion of the K -term does not sub- 
stantially modify the other model parameters derived by 
least squares fit. 

It is very difficult to find in the literature other esti- 
mations of K at these large heliocentric distances, since 
in the majority of cases authors consider an axisymmet- 
ric rotation curve. But some authors have pointed out the 
persistence of a residual in the radial velocity equations. 
Comeron & Torra ( |l994j ) found K = -1.9 ± 0.5 km s" 1 
kpc" 1 from 0-B5.5 stars and K = —1.3 ± 0.9 km s" 1 
kpc -1 from B6-A0 stars with R < 1.5 kpc. The radial ve- 
locity residual for Cepheids was first recognised by Stibbs 
( 1956| ). Pont et al. (1994) studied three possible origins for 
it: a statistical effect, an intrinsic effect in the measured 
radial velocities for Cepheids (7-velocities) and a real dy- 
namical effect. They finally suggested that it could be due 
to a non-axisymmetric motion pro duced by a central bar 
of « 5 kpc in extent. Metzger et al. ( |l998[ ) found a residual 
of « —3 km s" 1 for a sample of Cepheids when consid- 
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ering an axisymmctric galactic rotation. They concluded 
that it might be due to the influence of spiral structure, 
not included in their model. However, in the present work 
we found a non-null K value even taking into account 
spiral arm kinematics. An understanding of the physical 
explanation of the AT-term requires further study. 



Appendix A: systematic velocity components in 
the proposed galactic model 

In this appendix we show the expressions of the velocity 
components in the three systematic contributions consid- 
ered in our galactic model: solar motion, differential galac- 
tic rotation and spiral arm kinematics. 



8. Conclusions 



A. 1. Solar motion 



Hipparcos astrometric data were used to derive the spiral 
structure of the Galaxy in the solar neighbourhood. We 
considered two different samples of stars as tracers of this 
structure. In the sample of O and B stars the astromet- 
ric data were complemented with a careful compilation 
of radial velocities and Stromgren photometry, providing 
reliable distances and spatial velocities for these stars. In 
the sample of Cepheid stars we used two period-luminosity 



relations, one considering a short cosmic scale (Luri 2000) 
and another one standing for a large cosmic scale (Feast 
& Catchpole |i"997l ). 

A kinematic model of our galaxy that takes into ac- 
count solar motion, differential galactic rotation, spiral 
arm kinematics and a A"-term was adopted. The model 
parameters were derived via a classical weighted least 
squares fit. The robustness of our results was checked by 
means of careful insimulations. 

A galactic rotation curve with a value of A Oort con- 
stant of 13.7-13.8 km s _1 kpc -1 was found for the sam- 
ple of O and B stars. In the case of Cepheid stars, we 
found A = 14.9-16.9 km s _1 kpc -1 , depending on the 
case and the cosmic scale chosen. We confirmed the dis- 
crepancies appearing in A when samples with a different 
distance horizon were used (Oiling & Mcrrificld 1998). For 
both cosmic scales we found an acceptable coherence be- 
tween radial velocity, proper motion and combined solu- 
tions. Concerning the second-order term of the rotation 
curve, we always found a low value, compatible with zero. 

The study of the residuals for radial velocity data made 
it evident that a AT-term was needed, which was found to 
be K = —(1.4-3.2) km s -1 kpc -1 for O and B stars and 
K = -(0.8-1.2) km s -1 kpc" 1 for Cepheids. Although 
small, an apparent compression motion seems to exist in 
the galactic local neighbourhood for heliocentric distances 
up to 4 kpc, though its physical mechanism is still un- 
known. 

The phase of the spiral structure at the Sun's posi- 
tion obtained (?/>© = 284-20°) places it between the center 
and the outer edge of an arm. This is in good agreement 
with the spatial distribution of Cepheid stars. The angu- 
lar rotation velocity of the spiral structure was found to 
be ilp w 30 km s -1 kpc -1 , which places the Sun near the 
corotation circle. 
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acknowledges the FRD grant from the Universitat de Barcelona 
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A star with galactic longitude I and galactic latitude b has 
the following radial and tangential velocity components 
owing to solar proper motion: 

v Tl = — Uq cos I cos b — Vq sin I cos b — Wq sin b 
v [j = Uq sin I — Vq cos I 

«bi = Uq cos I sin b + Vq sin I sin b — Wq cos b (A.l) 

where Uq, Vq and Wq are the components of the solar 
motion in galactic coordinates. 

A. 2. Galactic rotation 

We consider axisymmetric differential rotation of our 
galaxy, with a rotation curve that can be developed in 
the solar neighbourhood as: 



@{w) « 9(n7 Q ) + 



ge 

dw 



Az 



d 2 e 

dw 2 



Q(wq) + a T Aw + 6 r Aa7 2 



(A.2) 



where Aw — w — Wq (wq is the galactocentric distance of 
the Sun) and Q(wq) is the circular velocity at the Sun's 
position. We note that a r allows us to calculate the A Oort 
constant in the Sun's vicinity: 



1 


~<d(w G ) 




1 


2 


Wq 


(<9w) Q _ 


~~ 2 



7 o 



(A.3) 



The radial and tangential velocity components of a star 
in the galactic plane due to differential galactic rotation 
are: 

u r2 = ©C^ 7 ©) [sin(Z + 9) — sin/] cos 6 
+a r Azv sm(l + 9) cos b 
+b r Aw 2 sin(l + 8) cos b 



v h = ®( w q) [cos(Z + 9) — cos/] 
+a T Aw cos(Z + 9) 
+6 r Atzj 2 cos(Z + 9) 

v\ i2 = —Q{wq) [sin(Z + 9) — sinZ] sin b 
— a r Aw sin(Z + 9) sin b 
-b x Aw 2 sin(Z + 9) sin b 

where 9 is the galactocentric longitude of the star. 



(A.4) 
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A. 3. Spiral structure kinematics 



Lin's t heory (Lin & Shu 1964 Lin et al. 1969 see also 



Rohlfs 1977 ) assumes a spiral potential of the form: 



Vb = Acostfj 



(A.5) 



(.4, is the amplitude - A < - and tp is the phase of 
the density wave) which disturbs the axially symmetric 
gravitational potential. The shape of the spiral arms is 
well represented by a logarithmic spiral: 



(A.6) 



The amplitude q is a slowly varying function of w and the 
phase of the spiral structure can be related to the phase 
at the Sun's position by the following expression: 



ro (t) 



i)(vj,9,t) = v©O)+ra(n p t-0) + 

tanz 

mln^- 

=>V(w,M = 0) = ip Q +m9 + ^- (A.7) 

tani 

where f2 p is the angular rotation velocity of the spiral pat- 
tern, m the number of spiral arms and i the pitch angle (for 
trailing spiral arms, i < 0). The phase of the spiral struc- 
ture at the Sun's position and the pitch angle can be de- 
termined from optical and radio indicators. Nevertheless, 
we point out that the maximum in the distribution of spi- 
ral arm tracers (position of the observed spiral arms) may 
be shifted in relation to the minimum in t he pe rturbation 
potential (defined as ip = 0°; see Roberts 1969 ). 

The mean peculiar velocities due to the spiral arm per- 
turbations on the velocity field are, in the gas approxima- 
tion, the following: 



IIi 
9i 



kA 



K 1 — V 2 

1 kAw 



■ cos i/j = Lib cos ip 



1 



2 e i 



■ sin-0 = —0b sin-0 



(A. 



IT is positive towards the galactic anti-center and 0i is 
positive towards the galactic rotation. In tightly wound 
spirals (i.e., those with | tani| <C 1) the amplitudes Lib and 
0b are slowly varying quantities with the galactocentric 
distance, k is the radial wave number (for trailing spiral 
arms, k < 0): 



k = _A_ I ™ 

dw \ tani 



mln-^- 



m 



vutani 



(A.9) 



v is the dimensionless rotation frequency of the spiral 
structure, expressed in terms of the epicyclic frequency 
(«): 



ml e 
v = — SL 



(A.10) 



(notice that v < in the region with f£ p < fl = Q/zu, i.e. 
inner to the corotation circle). Furthermore: 



29 2 



vj dQ 
O dvo 



(All) 



and x is the stability Toomre's number (Toomre 1969) 
defined as: 

x = ^ (A.12) 



where a is the velocity dispersion of the gas particles. 
Since the velocity amplitudes Lib and Ob depend on this 
velocity dispersion, we introduce a dimensionless param- 
eter (/q) that relates the velocity amplitudes of the Sun 
to those of the sample stars: 



n 



1 



b0 



^stars 



1 



EL 



e 



i 



bo 



+ Xq 



/ Q n b 



^e b = / e b 



(A.13) 



1 - V 2 + Xq 

The so-called Lindblad resonances are defined as: 

(A.14) 



where the — sign corresponds to the inner resonance and 
the + sign to the outer one. In the region between both 
resonances (|^| < 1), 8b is always positive and Lib bas a 
sign that depends on the sign of v. 

The amplitude of the spiral potential can be expressed 

as: 



A 



K n h l 



(A.15) 



The maximum value of the radial force owing to this spiral 
potential is: 

F™ ax = k\A\ (A. 16) 

whereas the radial force due to the axisymmetric field is: 
dV 



F ^ = ^~ K^q (A.17) 
Therefore, the ratio between both quantities is: 
F£** Kil b l-v 2 +x 



/r = 



•FrO 



(A.18) 



The contributions in the radial and tangential velocity 
components of a star due to the spiral arm perturbation 
velocities IIi and 9i arc the following: 

v l3 — — LIi cos(Z + 9) cos b + U 1& cos I cos b 
+6i sin(/ + 9) cos b — 6i sin I cos b 
= — LI b cosipQ 

( In- 



x cos 



"0 



tan i 



-Lib sin sin 

-6 b sin ipQ 



cos(Z + 9) — f & cos I j cos b 
cos(Z + 9) cos b 



ro 
tani 



x cos 



In -2-* 

ro 

tani 



sin(^ + 9) — / Q sin I J cos b 



+6b cos tpQ sin 



e- ro ° 



tan i 
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v\ 3 — 111 sin(Z + 6*) — IIiq sinZ + 
9i cos(Z + 6) — 0i0 cos Z 

= lib COS IpQ 

( In 
m\9- 



X ^cos 

+IIb sin ipQ sin 
-6 b sin il> Q 
x ^cos m ^9 - 

+6b cos tpQ sin 



ro 

tan z 



m \ V - 



sin(Z + 0) — /© sinZ 
sin(Z + 9) 



A3 


= w G 




= e(G7 ) 


as 


= a r 


a 6 


- b r 


a 7 


= lib COS tpQ 


a 8 


= nbSinV'© 


a 9 


= b sin^ 


aio 


= Gb cos -00 



In-=- 

ro 

tan i 



In- 



tan i 



cos(Z + 0) — /© cos Z 
cos(Z + 0) 



Ln-2- 

ro 

tani 



«b 3 = III cos(Z + 6*) sin 6 — IIiq cos Z sin 6 
—01 sin(Z + 9) sin b + 0i Q sin Z sin 6 
= lib cos V'q 



(A.21) 

and fj(R,l,b) are functions of the heliocentric distance 
and the galactic longitude and latitude: 

f[ = — cos Z cos b 

f 2 = — sin Z cos b 

fl = -sin b 

fl = [sin(Z + 9) — sin Z] cos 6 

/g = Atu sin(Z + 9) cos 6 

fl = Aw 2 sin(Z + 9) cos 6 



x cos 



m 



In -2- 

^0 

tan i 



+IIb sin ipQ sin 

+© b sin-!/'Q 



In 



cos(Z + 9 — /q cos Z J sin 6 
cos(Z + 9) sin 6 



fr 
fl 



cos 



sin 



ro 



tan i 



cos 



m 



In-=- 

ro0 

tani 



In -2- 

ro 

tani 



In -2- 

ro 

tan i 



cos(Z + 0) — /© cos Z J cos 6 
cos(Z + 9) cos 6 
sin(Z + 9) — /© sin Z cos 6 



x cos 



In- 



TO 



tan i 



sin(Z + 9) — /© sinZ sin 6 . r 
/ /io 



sin 



— ©b cos ipQ sin 



m 



tan z 



In- 



tan i 



sin(Z + 0) cos 6 



sin(Z + 6 1 ) sin b 

(A.19) 



/} = sinZ 

fl = -cosZ 
/s = 



A. 4. Systematic velocity field in the proposed galactic ^ cos ^ + ^ cos ^ 

f 5 = Awcos(l + 6) 



model 

The systematic radial and tangential velocity components 
of a star in our galactic model are given by: 



fl = Am 2 cos(l + 9) 
( In- 



Y.ajf^R^b) 

10 

ui = v h + v h + v h = ^2 a jfj( R > ^ b ) 



f 7 = cos 

fl = sin 



m V - 



10 



(A.20) 



/i = 



— cos 



«b = ^bi + «b 2 + «b 3 = ^2 a,jfj{R, I, b) 

3 = 1 

where the constants aj contain combinations of the kine- /{ = sin 
matic parameters that we wish to determine: 



tani 

In -2-' 

ro 

tani 
ki- 



rn w - 



OT0 



tan i 



In 

^0 
tan z 



sin(Z + 0) — /© sin Z 
sin(Z + 0) 

cos(Z + 6*) — /© cos Z 
cos(Z + 9) 



oi = Z7© 

«2 = V© 



fl = cos Z sin Z) 
/ 2 b 



/a = sin Z sin 6 
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f 3 b 



ft 



- cos b 

- [sin (7 + 6*) — sin I] sin 6 
-Aw sin(Z + 6) sin 6 
-Azn 2 sin(7 + 6) sin 6 



The distance error of the pseudo-star has a distribution 
law: 



e(R) 



(B.l) 



/ 8 b = 



Jo 



f b 

1 10 



cos 



tani 



tani 



tani 



cos(Z + 0) ~ f Q cos Z 1 sin 6 



cos(Z + 0) ; 



where <r# is the individual error in the photometric 
distance of the real star (Rq). 

To generate the kinematic parameters we randomly as- 
signed to each pseudo-star a velocity (U, V, W) by assum- 
ing a cosmic dispersion (au, cry, aw) and a Schwarzschild 
distribution: 

i( u-Uq \ 2 i( v-Vq \ 2 i( w-w \ 2 

ip' v (U, V,W) = e n "v ) »l "v J A "w ) 



In- 



tani 



In our resolution procedure the parameters a,j are com- 
puted followi ng an iterative scheme (extensively explained 
in Fernandez 199S ) and, from these, the kinematic param- 
eters Uq, Vq, W@, a T , b Y , ipo, lib, ©b and /© are derived. 

With regard to the free parameters of our model, dif- 
ferent values for m, i, m®, Q(wq) were considered (see 
Sect. |7h. 



Appendix B: simulations to check the kinematic 
analysis 

In Paper I we did numerical simulations in order to evalu- 
ate the biases in the kinematic model parameters (in that 
case, the Oort constants and the solar motion compo- 
nents) induced by our observational constraints and er- 
rors. In the present work, we have also generated sim- 
ulated samples in the same way, though the significant 
correlations detected between some parameters make it 
advisable, in this case, to carry out a more detailed study. 

In this section we present the process used to generate 
the simulated samples (the same as in Paper I, except for 
the change in the systematic contributions considered), 
the results we obtained and, finally, the quantification of 
the biases present in our real samples. 

B.l. Process used to generate the simulated samples 

To take into account the irregular spatial distribution of 
our stars and their observational errors, parameters de- 
scribing the position of each simulated pseudo-star were 
generated as follow: 

— From each real star we generated a pseudo-star that 
has the same nominal position (Ro,l,b) -not affected 
by errors- as the real one. 

— We assumed that the angular coordinates (l,b) have 
negligible observational errors. 



(B.2) 

sin(Z + 6) — / Q sin/ sinfc where (U , Vq, W ) are the reflex of solar motion. These 
/ components were transformed into radial velocities and 

proper motions in galactic coordinates using the nominal 
sin(Z + 6) smb (A.22)position of the pseudo-star (Rq, I, b). The systematic mo- 

tion due to galactic rotation and spiral arm kinematics 
was added following Eqs. (A. 4) and (A. 19), obtaining the 
components (v Io , jj,\ , jjb\, ) for each pseudo-star. Finally, in- 
dividual observational errors were introduced by using the 
error function: 



e(tt r ,Mi)M>) = e 



a,., and a,, 



(B.3) 

are the observational errors of the 



where a Vl 
real star. 

At the end of this process we had the following data 
for each pseudo-star: galactic coordinates (R,l,b), veloc- 
ity parameters (v T , Mi, Hb), errors in the velocity parame- 
ters (a Vi , er^j , ) and error in the photometric distance 
(o\r). The simulated radial component of those pseudo- 
stars generated from a real star without radial velocity 
was not used, thus we imposed on the simulated sample 
the same deficiency in radial velocity data that is present 
in our real sample (see Sect. 2.2 and Appendix B in Paper 
I for more details). 

Following this scheme, several sets of 50 simulated 
samples for both O and B stars and Cepheids were built. 
A classical solar motion of (U, V, W) = (9, 12, 7) km s -1 
was considered, taking the dispersion velocity compo- 
nents (ajj ,ay ,uw) = (8,8,5) km s _1 for O and B stars 
(see Paper I) a nd (a y, ay, aw) = (13,13,6) km s" 1 for 
Cepheids (Luri 2000| ). For the galactic rotation parame- 
ters, we chose the values a r = —2.1 km s" 1 kpc" 1 and 
b r = 0.0 km s _1 kpc -2 , which correspond to a linear rota- 
tion curve with an A Oort constant of 14.0 km s" 1 kpc -1 . 
On the other hand, for the spiral structure parameters 
several sets of values were used for ip G (from tpQ = 0° to 



315°, in steps of 45°) and f2 p (from Sl p = 10 km 



s" 1 kpc -1 to n p = 40 km s" 1 kpc" 1 , in steps of 5 km s" 
kpc" 1 ), whereas a fixed value of / r = 0.05 was considered 
(Yuan 1969). From (au, cry, aw), Or, and / r , the values 
of nb, Ob and / were inferred for each set of samples. 

56 sets of 50 samples for both O and B stars and 
Cepheids were generated. Concerning the free parameters 
of our model, in a first stage we adop ted classical val- 

1969| wq — 8.5 kpc, 



ues 



(rn = 2, i = —6°, Lin et al. 
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9 km s~ x 




12 km s- 1 


Wq 


7 km s" 1 


(<tu , <?v , crw) 


(8, 8, 5) km s^ 1 (O and B stars) 




(13, 13, 6) km s~ (Cepheid stars) 


Hi 


-2.1 km s" 1 kpc" 1 




0.0 km s" 1 kpc" 2 



fr 0.05 

ip @ from 0° to 360°, 

in steps of 45° 

Q p from 10 km s _1 kpc -1 to 40 km s -1 kpc -1 , 

in steps of 5 km s kpc 



9(ro Q ) = 220 km s" 1 , Kerr & Lyndell-Bcll |1986[ ), though 
we also tested cases with m = 4, i = —12° (Amaral & 
Lepine |1997 ) and Wq = 7.1 kpc, Q(zuq) = 184 km s -1 



(Oiling & Merrifield 1998). In Table B.l we summarize all 
the adopted kinematic parameters. 

B.2. Results and discussion 

B.2.1. Results for a 2-armed model of the Galaxy 

A complete solution simultaneously taking into account 
radial velocity and proper motion data was computed. Our 
test showed that the number of Cepheids within 2 kpc 
from the Sun prevents the obtainment of reliable results. 
In Figs. p3.l| and B.2 we show the results obtained for the 
simulated samples of O and B stars (0.6 < R < 2 k pc) and 
Cepheids (0.6 < R < 4 kpc) in case A (see Table |Bj]). 

As a first conclusion, and confirming our suspicions, a 
systematic trend with tpQ and/or Q p is observed in most 
cases. This behaviour is produced by the correlations be- 
tween some terms in the least squares fit, which depends 
on the spatial distribution of each sample. 

For solar motion a bias between —1.5 and 1.5 km s _1 
(depending on ipQ and Q p ) was found for Uq and Vq, and 
of only —0.3 km s _1 for Wq. For both O and B stars and 
Cepheids we found the bias on Vq and Wq to be indepen- 
dent of fi p , with a slight dependence on ipQ . For i/;q = 270° 
and fip = 10 km s _1 kpc -1 we found a large negative bias 
on Vq , but with a great standard deviation. This occurred 
in several samples (inside this set) with serious conver- 
gence problems in the iteration procedure we use to solve 
the least squares fit. Similar problems in other cases with 
V>© = 270° will be found later. For O and B stars, the 
standard deviations in the solar motion components are 
w 0.6 km s _1 for Uq and Vq (except for ipQ — 270°), and 
« 0.3 km s _1 for Wq. On the other hand, in the case of 
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Fig. B.l. Bias (obtained value— simulated value; left) and 
standard deviation (right) for each set of 50 samples in the 
solar motion (top), galactic rotation (middle) and spiral 
arm kinematic (bottom) parameters for O and B pseudo- 
stars (Case A). Values of fi p : 10 km s _1 kpc -1 (black 
double solid line), 15 (black solid line), 20 (grey solid line), 
25 (dotted line), 30 (dashed line), 35 (long dashed line) 
and 40 (dot-dashed line). 
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Fig. B.2. Bias (obtained value— simulated value; left) and 
standard deviation (right) for each set of 50 samples in the 
solar motion (top), galactic rotation (middle) and spiral 
arm kinematic (bottom) parameters for Cepheid pseudo- 
stars (Case A). Values of fi p : 10 km s -1 kpc -1 (black 
double solid line), 15 (black solid line), 20 (grey solid line), 
25 (dotted line), 30 (dashed line), 35 (long dashed line) 
and 40 (dot-dashed line). 



Cepheids these values increased to w 1.4-2.0 km s 1 and 
w 0.8 km s _1 , respectively. 

The biases found in the first- and second-order terms 
of the galactic rotation curve are negligible for Cepheids, 
with a level fluctuation of ±0.3 km s _1 kpc -1 (or km s _1 
kpc -2 ). In this case there is a standard deviation of 1.3 
for a r and 0.5 for b r . For O and B stars the biases clearly 
depend on varying from —0.7 to —0.5 km s" 1 kpc -1 
for a T , and from —1.0 to 0.1 km s -1 kpc -2 for b r . The 
standard deviations are 0.8 km s -1 kpc -1 and 1.2 km s -1 
kpc" 2 , respectively. 

Let us study the biases that have an effect on the de- 
termination of the spiral structure parameters. As a gen- 
eral conclusion, Figs. [B.l| and show that our sample 
of O and B stars supplies better results than the Cepheid 
sample. 

In the case of O and B stars, we found a clear depen- 
dence with ?/>0 and tt p in ip G , lib and 6b determinations, 
whereas /© and fl p only show peculiar behaviour around 
ipo = 270°. Concerning tpQ, the bias oscillates from —20° 
to 30°. The standard deviation of the mean for the 50 
samples of each set is about 10-20°. On the other hand, 
for the amplitudes lib and Ob the biases are of ±2 km 
S _1 , with a standard deviation of about 1 km s _1 . Neither 
/© nor Hp have a considerable bias, except for ?/>© = 270°, 
where both biases and standard deviations go up. 

For Cepheid stars similar results were obtained, but 
with larger standard deviations in all cases. The bias in 
■0© changes from —20° to 10°, with standard deviations 
of about 30-60°. In the case of lib, we found a clear de- 
pendence on both -00 and fl p , with a bias of ±1.5 km s -1 
and a standard deviation of about 2 km s _1 . On the other 
hand, for Ob the bias is smaller, from to 0.8 km s _1 , 
and the standard deviation is 1-1.5 km s~ x . As for O and 
B stars, for /© and f2 p small biases are found, though the 
standard deviations are larger in this case. 



B.2.2. Results considering possible errors in the choice 
of the free parameters 

An interesting point to analyse is the study of the biases 
produced by a bad choice of the free parameters in our 
model (m, i, vjq, 0(c*7©)). In the same way as in the 
previous section, we simulated 50 samples for each one of 
the cases considered in t he real resolution, i.e. cases A, B, 
C and D (see Table B.2). The simulated parameters were 
the same as in Table B.l for solar motion and galactic 
rotation. For spiral arm kinematics, we considered ipQ = 
315° and il p — 30 km s _1 kpc -1 (similar values to those 



obtained from real samples; see Sect. ||). 



In Tabic B.2 wc show the biases and standard devia- 
tions when solving the model equations in crossed solu- 
tions (e.g. we generated 50 simulated samples considering 
the free parameters in case A, and then we solved equa- 
tions using the free parameters adopted for cases A, B, C 
and D, and so on for the other cases). As a first conclusion, 
we can observe that a bad choice in the free parameters 
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Table B.2. Bias and standard deviation in ip Q and fl p for crossed solutions for the simulated samples of and B 
stars and Cepheids. Units: tp@ in degrees; Q p in km s _1 kpc" 1 . 





and B stars with 0.6 < R < 2 kpc 


Cepheid stars with 0.6 < R < 


4 kpc 




Case A 


Case B 


Case C 


Case D 


Case A 


Case B 


Case C 


Case D 


Case A simulated 




23. 


17. 


24. 


19. 


-11. 


0. 


-2. 


3. 




28. 


27. 


25. 


25. 


65. 


61. 


71. 


71. 




2.2 


3.1 


-1.1 


-0.8 


-4.3 


-3.5 


-2.5 


-4.0 




6.2 


7.3 


2.9 


3.5 


18.9 


16.6 


8.8 


7.3 


Case B simulated 


AV© 


32. 


24. 


32. 


25. 


0. 


4. 


-11. 


1. 




34. 


32. 


31. 


30. 


83. 


75. 


82. 


78. 


AQ P 


0.7 


1.9 


-1.6 


-1.3 


-4.3 


-14.3 


-4.3 


-4.1 




7.1 


8.0 


3.8 


3.9 


16.2 


52.8 


7.0 


8.1 


Case C simulated 


A^© 


18. 


14. 


17. 


14. 


4. 


2. 


-2. 


10. 


a i>e 


24. 


23. 


21. 


21. 


62. 


58. 


60. 


63. 




6.7 


8.2 


1.2 


2.0 


-2.4 


-1.1 


-1.0 


-3.0 




6.5 


7.3 


2.9 


3.3 


25.1 


21.5 


10.9 


9.2 


Case D simulated 


AV>© 


29. 


23. 


25. 


21. 


9. 


10. 


-9. 


1. 




31. 


29. 


27. 


26. 


84. 


78. 


81. 


73. 


AQ P 


4.3 


6.3 


0.2 


1.1 


-3.9 


-4.9 


-3.2 


-2.8 




7.1 


8.2 


3.7 


3.8 


17.1 


19.5 


7.8 


8.6 



does not substantially alter the derived kinematic param- 
eters, particularly ipQ. In other words, for each set of sim- 
ulated samples we obtained nearly the same values for the 
parameters whether we solved the Eqs. (|J) with the cor- 
rect set of free parameters or with a wrong combination 
of them. Differences in -0 Q do not exceed 10° for O and 
B stars and 20° for Cepheids. In the case of tt p we found 
large differences in some cases, but always when the stan- 
dard deviation was also large. This is especially true for 
Cepheids. A remarkable point is that the minimum bias 
was not always produced when we properly chose the free 
parameters. 



B.2.3. Conclusions 



In the light of these results, we conclude that we are able 
to determine the kinematic parameters of the proposed 
model of the Galaxy from the real star samples described 
in Sect. ^, supposing that the velocity field of the stars is 
well described by this model. We studied case A (m = 2, 
i = —6°, Wq — 8.5 kpc, 8(zn ) = 220 km s _1 in detail in 
these simulations, but we also looked at the other combi- 
nations of the free parameters (cases B, C and D), with 
similar conclusions. Nevertheless, the study of crossed so- 
lutions has shown that it will be very difficult to decide 
between the several s et o f free parameters discussed in 
Sect. ^ (see also Table jU2), owing to the small differences 
obtained when changing the free parameters in the condi- 
tion equations. 
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